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I present results for the two-loop self-energy functions for neutral and charged Higgs scalar bosons 
in minimal supersymmetry. The contributions given here include all terms involving the QCD cou- 
pling, and those following from Feynman diagrams involving Yukawa couplings and scalar interac- 
tions that do not vanish as the electroweak gauge couplings are turned off. The impact of these 
contributions on the computation of pole masses of the neutral and charged Higgs scalar bosons is 
studied in a few examples. 
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I. INTRODUCTION 

The small ratio of the electroweak symmetry breaking 
scale to other possible energy scales, including the Planck 
scale, is one of the most important puzzles in high-energy 
physics today. This hierarchy is stabilized by low-energy 
supersymmetry [1-3], but only if it is within discovery 
reach of the Large Hadron Collider (LHC) and subject to 
detailed study at a future TeV-scale e + e~ linear collider 
(LC). It is a pleasant feature of supersymmetry that the 
Higgs sector is both pcrturbatively calculable, and highly 
sensitive to radiative corrections at least at two-loop or- 
der. The precision of measurements at the next gener- 
ation of high-energy physics experiments will therefore 
allow precision tests of theoretical model frameworks for 
low-energy supersymmetry. For example, the mass of the 
lightest neutral Higgs scalar boson, h°, may be obtained 



at the LHC with an uncertainty of perhaps 100-200 MeV 
[4], and about 50 MeV at a LC [5]- [7]. 

Much work (see for example [8]- [37] and references 
therein) has already been done on radiative corrections 
to electroweak symmetry breaking in supersymmetry and 
on the related problem of evaluating the physical masses 
of the Higgs scalar bosons. The effective potential for 
the minimal supcrsymmetric standard model (MSSM) 
has now been evaluated at two-loop order [30, 31]. In 
addition to allowing an accurate implementation of elec- 
troweak symmetry breaking in the MSSM, this allows a 
calculation [32] of the physical mass of the lightest neu- 
tral Higgs scalar, incorporating the complete one-loop 
results and all two-loop results in the effective potential 
approximation. This means that the two-loop self-energy 
contributions to the pole mass are estimated by setting 
the external momentum invariant equal to zero, instead 
of evaluating them at the pole in the rcnormalizcd propa- 
gator. This can be a good approximation for the lightest 
Higgs scalar, since it is much lighter than most of 
the virtual particles propagating in loops, in particular 
the top quark and the squarks. However, the error made 
in doing so is still significant compared to the eventual 
experimental uncertainty in the mass of a light Higgs 
scalar boson, as we will see below. Also, it is gener- 
ally not a valid approximation for the other Higgs scalar 
bosons if , A , fi^, especially in the limit that they are 
heavy. In order to adequately compete with the accu- 
racy that can be obtained at the LHC and an LC, it will 
probably be necessary to have the complete momentum- 
dependent set of corrections to the two-loop self energy, 
and the leading three-loop contributions in the effective 
potential approximation, at least for h . 

In this paper, I extend previous work by presenting an- 
alytical expressions for some of the leading contributions 
to the two-loop self-energy functions for the Higgs scalars 
in the minimal supcrsymmetric standard model. My 
calculations use the mass- independent DR scheme [38] 
based on regularization by dimensional reduction [39] . Of 
course, these results should eventually agree with calcu- 
lations done in the on-shell type schemes for all questions 
posed in terms of physical observables, up to corrections 
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of higher order. As a matter of opinion, I find the calcula- 
tions in the DR scheme to be simpler than in the on-shell 
schemes, and more flexible in the sense that they can be 
performed once for generic field theories and then ap- 
plied to all kinds of special cases. (In on-shell schemes, 
the organization of the calculations depends on a special 
choice of observable input parameters; this choice will 
be different for different particles and for different the- 
ories.) Indeed, the results presented below rely on cal- 
culations already performed in a generic renormalizable 
quantum field theory in ref. [40]. In that paper, formu- 
las for the two-loop scalar self-energy diagrams involving 
up to two gauge couplings were presented in terms of a 
minimal basis of two-loop integrals. Explicit definitions 
and procedures for the efficient numerical evaluation of 
these basis integrals 1 are described in ref. [41, 42]. Com- 
parisons with the predictions of specific models for very 
high-energy physics and supcrsymmctry breaking will re- 
quire the evaluation of DR scheme parameters, by global 
fits of many observables to data. 

The objects of interest in this paper are the one-loop 
and two-loop contributions to the self-energy function 
matrices for Higgs scalar fields cpf. 
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as functions of the squared-momentum invariant 
s = —p 2 . 



(1.2) 



using a metric of signature ( — h++). Here s is always 
given an infinitesimal positive imaginary part to resolve 
branch cuts above thresholds. Then the gauge-invariant 
[48]-[51] complex pole masses of the Higgs scalar bosons, 



s k = Mi - iT k Mk, 



(1.3) 



can be found by iteratively solving the equation 

Dot [(to 2 - s^Sij + Uij(a k )] = 0, (1.4) 

where the mf are the tree-level renormalized running 
squared masses. Here, the self-energy function must be 
evaluated in the sense of a Taylor series around a nearby 
point on the real s axis; in other words, the self-energy 
and its derivatives are first evaluated for s with an in- 
finitesimal positive imaginary part, and this data is then 
used to construct a Taylor series expansion for complex 
s. This is necessary because the imaginary part of the 
pole mass is negative, while the imaginary part of s is al- 
ways positive. One representation of the solution, which 



1 The basis integrals are renormalized versions of the ones whose 
recursion relations were worked out in [43] and implemented in 
[44]. The strategy for their evaluation in [41] (soon to be im- 
plemented in a computer program package [42]) is similar to the 
one put forward earlier in [45]. Some other useful two-loop self- 
energy basis integral strategies are found in [46]- [47]. 



maintains manifest gauge invariance at each order in per- 
turbation theory, is 



Det 



(to 2 - s k )6ij + [U k ]i 



0. 



where, at one-loop order, the solution is obtained 
using 



(1.5) 
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and then at two-loop order, 
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[Uk]ij = i^^VD + T^^-^n^VD 
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(16tt 2 )< 



ngV 2 ). 



(1.7) 



Formally, the difference between this method and the 
method of iterating eq. (1.4) directly is of three-loop 
order. However, the tree-level value of to 2 ^ runs quite 
rapidly with the rcnormalization scale Q, so perform- 
ing a Taylor series expansion about it is formally valid 
but numerically suspect. The difference between these 
two methods for computing the pole masses of the Higgs 
scalars usually turns out to be small for the real parts, 
but the procedure of iterating eq. (1.4) directly gives a 
result for the imaginary part of the complex pole mass of 
h° that is much more stable with respect to changes in 
the renormalization scale Q. 

The calculations used in this paper neglect the Yukawa 
couplings of the first two families, and the corresponding 
soft (scalar) 3 interactions. Thus, I use as inputs the fol- 
lowing 33 DR parameters at a specified renormalization 
scale Q: 



V u , Vd, 

93, g, g', vu Vb, y T , 

m Q t , m u, m 2 u ., to 2 ,., to 2 



e« ' 



(2 = 1,2,3) 



l H d , 



b, n, 



M 3 , M 2 , Mi, a t , a b , a T , 



(1.8) 
(1.9) 
(1.10) 
(1.11) 
(1.12) 



in the notation of refs. [3, 31]. No assumptions regarding 
CP-violating phases are made, so the last 7 parameters 
may be complex. The other parameters are always real, 
either by definition or by convention, without loss of gen- 
erality. This means that the formulas below are valid for 
general CP violation in the soft terms of the MSSM, but 
neglecting the usual Cabibbo-Kobayashi-Maskawa CP vi- 
olating parameter. At tree-level, there is no CP-violation 
in the Higgs sector, so one defines tree-level mass eigen- 
states 0° = (h°,H ,G°,A ) and 4>f = {G ± ,H ± ) with 
the usual CP quantum number assignments. In general, 
the self-energy functions then consist of a 4 x 4 matrix 
for the neutral scalars and a 2 x 2 matrix for (j>f. 
The parameters v u and Vd are actually redundant; they 
are defined to be the Landau gauge vacuum expectation 
values of the Higgs fields at the minimum of the two-loop 
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effective potential evaluated at Q. In practice, they can 
be taken as given and used to eliminate b and or vice 
versa. 

In calculating the effective potential, and the self- 
energies below, I use the Landau gauge for electroweak 
bosons, and a general covariant gauge for gluon propa- 
gators. The fact that v u and v d minimize the Landau 



gauge two-loop effective potential means that the sum 
of all tadpole diagrams, including the tree- level contri- 
butions, vanishes identically through the same order, so 
that they do not need to be included explicitly in pertur- 
bative calculations. As in ref. [31], the tree-level neutral 
Higgs squared mass matrices are therefore given by: 



f\^\ 2 + m 2 Hu + (9 2 +9 ,2 )(3vl^vj)/4 -b (g 2 + g' 2 )v u v d /2 \ 

m K ~ \ -b- (g 2 + g' 2 )v u v d /2 | M | 2 + m\ d + (g 2 + g' 2 )(3v 2 v 2 u )/a) ' 

2 _ f\^\ 2 + m 2 Hu + (g 2 + g' 2 )(v 2 u -v 2 d )/4 b \ 

m V - { b W+m 2 Hd + {g 2 +g' 2 )(vl-vl)/A)' ^ U > 

in the (Re[ff°], Re[H d }) and (Im[iJ"], Im[_ff°]) bases, respectively. Note that even in the presence of arbitrary CP 
violation, the tree-level squared mass matrices always seperate into 2x2 blocks in this way, because of the freedom to 
choose b real and positive at any given renormalization scale. The complex charge ±1 Higgs scalar tree-level squared 
masses are obtained by diagonalizing the matrix 

(W + m\ + (g 2 + g' 2 )vl/A + (.g 2 - .g'> 2 /4 b + g*v u v d /2 \ 

^ ~ \ b + g*v u v d /2 H 2 +^ d + (. 9 2 + 3 '> 2 /4+( 3 2 - 3 '> 2 /4j' 



in the (H+,H d *) basis. 

There is another approach (see for example Appendix 
E of ref. [20]) in which the condition of vanishing of the 
tadpoles is used to replace the tree-level masses with 
different expressions, by eliminating |/i| 2 + m 2 H and 
|/i| 2 + m H in favor of combinations of b and m z . In 
that approach, tadpole terms do appear explicitly in the 
loop-level part of the mass matrices. Of course, both ap- 
proaches must agree in principle on their predictions for 
the physical masses, through whatever loop order one 
is working. In the approach followed here, the tree- 
level eigenvalues of the lightest neutral Higgs and the 
Goldstone bosons as obtained from eqs. (1.13)-(1.15) are 
rather strongly dependent on the choice of renormaliza- 
tion scale. These tree-level masses enter into the kine- 
matic loop integral functions. However, as we will see 
below, the resulting scale dependences of the calculated 
physical Higgs scalar masses are very small. A wide range 
of renormalization scales gives consistent results for the 
physical masses, within the uncertainties inherent in the 
two-loop approximation. (It is also possible to expand 
the analytical formulas presented here around any choice 
of tree-level squared masses, treating the differences as 
perturbations. The one-loop integral functions are all 
known analytically, so this does not present any techni- 
cal difficulties, but will not be explored in detail here.) 

In this paper, I include all one-loop corrections to the 
Higgs scalar boson self-energies. The two-loop correc- 
tions that are included are of two types. First, I include 
all diagrams that involve the QCD coupling g 3 . This 



includes all effects of order: 

22 2 22 22 2 / 2 /2 1c i 

9zVti g^ytVb, g 3 y b , g 3 g , g 3 gg , 333 > l 1 - 16 ) 

and those related by replacing one or both powers of y t 
or yb by the corresponding soft coupling at or This 
means all diagrams involving the gluon or the gluino, and 
also the diagrams involving the four-squark interactions 
proportional to g\. Second, I include all diagrams that 
do not vanish when the electroweak gauge couplings are 
turned off. These include effects proportional to 

yt, VtVb, yhh vtvl, yt yt, vlvl-, (i-i?) 

and those related by replacing one or more Yukawa cou- 
pling^) by the corresponding soft terms at,ab,a T . Also, 
I include electroweak effects whenever they contribute to 
the same Feynman diagrams as just mentioned. This in- 
cludes both explicit factors of g, g' in the scalar couplings 
that also involve Yukawa couplings, and implicit factors 
in the mixing angles of the Higgs scalars, squarks, slep- 
tons, neutralinos and charginos. (It would seem counter- 
productive to try to disentangle the latter anyway.) In 
the future when all of the two-loop self-energy contribu- 
tions become available, it will just be a matter of adding 
in the contributions of the Feynman diagrams not consid- 
ered here. It follows that some, but not all, effects of or- 
der e.g. y\g 2 are included in the present paper. Thus, the 
formal level of approximation is to neglect electroweak ef- 
fects not involving g 3 at two-loop order; but for future 
convenience some of them are included anyway. In the 
case of the lightest Higgs scalar boson h°, all other two- 
loop corrections to the self-energy are included in the 
effective potential approximation, as in refs. [30-32]. 
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In much of the parameter space of the MSSM, includ- 
ing the decoupling limit for the heavier Higgs scalars, the 
scalar h° is predominantly made out of the gauge eigen- 
state field that couples to the top quark, while H°, A , 
and are predominantly made out of the gauge eigen- 
state field that has a Yukawa couplings to the bottom 
quark. Therefore, because of the large top mass com- 
pared to the other quarks and leptons, the effects detailed 
above are generally more significant for h° than for the 
other Higgs scalars, at least when tan/3 is moderate. 

The Feynman diagram topologies that play a role in 
this paper are shown in Figure 1. Each diagram shown 
that involves fermions actually refers to several distinct 
ones, with chirality-reversing fermion mass insertions in- 
serted in all possible ways. For each diagram, there is a 
corresponding finite loop integral function, which also in- 
cludes the DR counterterms for that diagram. The label 
on each diagram refers to that function, strictly follow- 
ing the notation and definitions found in [40] , which lists 
them in terms of the minimal set of basis functions. The 
numerical evaluation of the basis functions is in turn de- 
scribed in detail in ref. [41]. 

The rest of this paper is organized as follows. Section 
II presents the complete list of three- and four-particle 
couplings used in the calculations. The known one-loop 
results for the Higgs scalar self-energies are reviewed in 
section III. The two-loop self-energy contributions de- 
scribed above are given for the neutral Higgs scalars in 
IV, and for the charged Higgs scalars in V. Section 
VI briefly recounts some consistency checks, and stud- 
ies some numerical results for specific model parameters. 

II. COUPLINGS 

In this section, I provide the list of three- and four- 
particle couplings needed in the rest of the paper. The 
conventions and notations for the MSSM Lagrangian pa- 
rameters and mixing matrices strictly follow those given 
in section II of [31], which will not be repeated here for 
brevity. 

[The signs of some of the couplings listed here do differ 
from those listed in section III of [31], namely equations 
(3.6)-(3.9), (3.11)-(3.13), (3.27)-(3.33), (3.35)-(3.38), and 
(3.44)-(3.47) of that paper. These sign conventions actu- 
ally make no difference at all for ref. [31], because three- 
particle couplings always appear squared in the two-loop 
effective potential. However, the signs are important in 
the present paper, and have been chosen to agree con- 
sistently with ref. [40]. To avoid confusion, all of the 
relevant couplings are listed here.] 

The couplings of fermions to the Higgs scalar bosons 
0° = (h a ,H°, G a ,A°) and <f>f = (G ± ,H ± ) are: 

Y ii<p° = ytk u <p°/V2, (2-1) 
Y bHl = Vbk d4> o/V2, (2.2) 
Y t-^ = Vrk d4> o/V2, (2.3) 



Y t b4 >t = -vtK^i ( 2 - 4 ) 

Y bt4,- = -yb k d4,+ , ( 2 - 5 ) 

Y t,J- = -Vrk d ^ t . (2.6) 
The fermion-neutralino-sfermion couplings are: 



Y - 

uNiU* L 


= 


(gN* 2 + 5 'V*/3)/V2, 


(2.7) 


Y - 

uN iUn, 


= 


-2V2 3 'V*/3, 


(2.8) 




= 


(-«+«/3)/\/2, 


(2.9) 


Y-~ ~ 

dNidji 




V2g'N* /3, 

V iJ III ' 


(2.10) 


Y -, 

± eN i e* L 




-(qN* 9 + q'N*-,)/V2, 


(2.11) 


Y - 






(2.12) 


Y - 




(gN* 2 - g'N*, )/V2, 


(2.13) 


Y - - 




Lt Y uN^ L +RtNkyu 


(2.14) 


Y-- - 

1 tN i t j 




Rt j Y m R + L h N ^Vu 


(2.15) 


Y bNib* 




L l Y dNd* L + R * b] N ayb, 


(2.16) 


Y-- - 

I bN i b j 




Rb Y dNd R + L b 3 N ^ 


(2.17) 








(2.18) 


Y - 






(2.19) 


lion-chargino-sfermion couplings are: 




Y - 

dCiU* L 




Y eCiO' - Y rC t i>; ~ 9 V ii, 


(9 90) 


Y - - 

udd* L 






(2.21) 


Y bC^ 






(2.22) 


Y-~ ~ 

I bC i t j 




- L iUtiVb, 


(2.23) 


Y tc,i; 




L%9Ua-RlUt 2 Vb, 


(2.24) 


Y-~ - 
I tC i b 3 




-L~ b V* 2 y u 


(2.25) 






L% j gXJ* x -R%U* 2 y T , 


(2.26) 


Y tCi9 t 




~VrU* 2 . 


(2.27) 



The neutralino and chargino couplings to Higgs scalar 
bosons are: 

Y c+c^l = + ( 2 - 28 ) 
Y N t N A o = (gN^-g'N^k^N*, 

-^°^)/ 2 + (*~i), (2.29) 
Y c t N A - k - K,t^U 

+V* 2 (gN* 2 +g'N* 1 )/V2], (2.30) 

Y c-n^ = *» t WiNU 

-U^gN^+g'N^/Vl}. (2.31) 

The couplings of electroweak gauge bosons to each other 
and to the Higgs scalar bosons are: 

9w^4>t = ^d^d^i - K^K<t>+ ) A (2-32) 
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FIG. 1: The one-loop and two-loop Feynman diagram topologies in this paper, by order of first appearance. Dashed lines are 
for scalars, solid lines for fermions, wavy lines for electroweak vector bosons and ghosts, and curly lines for gluons. The label 
on each diagram refers to a corresponding renormalized integral function, as defined in ref. [40]. There are 7 one-loop and 
40 two-loop distinct topologies here, accounting for fermion mass insertions (indicated below with a bar over the appropriate 
subscript F), but not counting separately diagrams obtained by exchanging external lines or reversing all fermion chiralities. 



9z4>°4>° = V9 2 +9' 2 I™[K<j,°K<t>° ~ fc d^^o]/2, (2.33) 
= l5 ^ 2 - .9' 2 )/ 2 vV+5' 2 , (2-34) 



9ww<j>fcf>j ~ ^io9 /2j 
9zz^° = %(.9 2 + .9 /2 )/2, 

5^ +0 7 =M 5 2 -g' 2 ) 2 /2(.9 2 + g' 2 ), 



(2.35) 
(2.36) 
(2.37) 
(2.38) 



9zz<f>° 

= (g 2 + g' 2 )Re[v u k u<t> o + v d k d ^o}/V2, (2.39) 

= 3 2 R-eK&„0° + v d k d4> o]/V2, (2-40) 

= e 9(vuk U(f ,+ - v d k d(j> +)/V2, (2.41) 

9wz<j>+ = - e 9'(v u k u<j> + - v d k d(f> +)/V2, (2.42) 

where e = g'g/ \/g 2 + g' 2 is the QED coupling. 



The couplings of four Higgs scalar bosons are given by: 

A ^^°^e = 2 + 9' 2 )R4K<t>°K<t,° ~ k d<p°k d ^} Re [ k u<l>° k k u^ m ~ k d<j> « k* d4>0 J/A + (i <-► k) + (i <-» m), 
= g 2 (SijS km + 2fc u0 ofc d0 ofc d0 +fc M0 + + 2fc* ofc* o/c u0 +fc d0 + ) 

+9' 2 {k u ^K^ ~ k d^k* d!t> o){k uit> +k uij)t - k dtf> +k d4> + )] /8 + (i <-» j), 



(2.43) 



A., 



?' 2 ' 

„2 , /2\ 



= (9 + 9 )(2fc *+k <+k , + k ,4 



^tt0+ 



k d<j>+ k ucj,+ k u4>+ 



(2.44) 

*«W+)/4+(«~d)- (2-45) 



and the couplings of three Higgs scalars are: 

^cfityi = (.9 2 + ff' 2 )Re[fc u0 ofc* o - k M ok* d4>0 ]Re[k u4 ,o k v u - fc d0 ou d ]/2\/2 + (k <-> i) + (fc <-> j) 
A 0>+0- = {5 2 (bdfc u< ^o + v u k drj> o]k d<f> +k u4i + + [v d k* U(f> o + v u k* d<j> o]k u<t> +k d<j} + + 6 jk Re[v d k d(j} o + v u k ulf> o}) 

+.9' 2 [^ +^+ - Krt k u<t>i]Mvdk d4 ,o - v u k u<j> o]}/2V2. (2.46) 

The couplings involving sfermions are conveniently written using the quantities Ij and Yj, defined to be the third 
component of weak isospin and the weak hyperchargc of the left-handed chiral superfield containing the squark or 
slepton /: 





u L 


d L 




eL 


UR 


d~n 






1/2 


-1/2 


1/2 


-1/2 











Y f 


1/6 


1/6 


-1/2 


-1/2 


-2/3 


1/3 


1 



Then we have for the couplings of two neutral Higgs scalars to sfermion-antisfcrmion pairs 

X <i>°<t>°ff* = (dfg 2 -Yfg /2 )Rc[k d4> ok d4> o - k u4 ,ok* U(j> o}/2 
for the sfermions / of the first and second families and v T , and 



(2.47) 



(2.48) 
(2.49) 
(2.50) 



Re[k d4> ok d(j) o]y b Skm + L b k L t m X 4>°i^d L di + R b k R t m X 4>°<j>°d R d H ' 

jTk?^ = Re [ k d^k d(j> o]y 2 5k m + L ?k L f m ^4>°4> ( >£ L e' L + R fk R f m ^<t>]e R e R 

for the other sfermions of the third family. The neutral Higgs-sfermion-antisfcrmion couplings are similarly given by 

A o//. = {Ifg 2 - Y f g' 2 )Re[k d4> av d - k u ^v u }/s/2 (2.51) 
for the first two families and v Tl and by 



X <t>°Jkt' m = ^v u y 2 Re[k u<p o}S km + L ik R* im (k ucf> oa t - k* d(j> o^*y t )/ V2 + Rf k L* im (K<j>» a t " k d^Vt)l V2 



+L t k L *i m X <l> iULui 



R t k R L X 4> ( ' 



'(U R U R ! 



(2.52) 



X 4>Tb k b- m = V2v d y 2 Re[k d4> o}S km + L^R^ {k d(f> oa b - fc^o/z* y b ) / 'V2 + R~ bk L~ m (k* d4> „a* b - k u4> o fiy b ) / V2 
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+L b k L i m x <f>°d L d* L + R b k R i m x <t>°d R d* R > ( 2 - 53 ) 

X 4>°i^f^ = V2v d y*Re[k d4 ,o]5 km + L fk R* fm {k dlt> oa T - k* u4>0 ^*y T )/V2 + R fk L* fm (k* d(t>0 a* - k ulf> ony T )/V2 

+ L f k L *f m X 4>°e L e* L + R fk R f m X 4>°e R e' R ( 2 - 54 ) 

for the other third family sfermions. The couplings of pairs of charged Higgs scalars to sfermions of the first two 
families are 

\Ujii- = (V + Y f9' 2 )(K,t k u,t - k d ,t k d,f)/ 2 - ( 2 - 55 ) 
For the sfermions of the third family, 

= R h R iM k u4>t k u^ + \t<t>-u R u R ) + L i k L tJy 2 b k dot k d4 + \uj™-r)> ( 2 - 56 ) 

ht^b k i* m = L b k L lJytK4>t k u^ + \t^-d L di) + R b k R i m (yl k d^ k d4>+ + ht^dnd-J' ( 2 - 57 ) 

= (y 2k d4>+ k d4>+ + V^*)' ( 2 - 58 ) 

V*7^r- = L r k L% m \ U -^ el + Rf k R*fJyr k d^ k d4>t + \Uj ^ ^ (2 ' 59) 

The charged Higgs-sfermion-antisfermion couplings are 

htdLul = = 9 2 (K^ + k d^d)/2 (2.60) 

for the first two families, and 

= L b k L L( x <j>+dLui - yt v uk u4> t - yl y dk d ^) - Ri k R *t m ytyb{k d<t> +v u + k u4> +v d ) 

- L b k R *t m (K<i,+ a t + kdt+^Vt) ~ R b k L *t m (k d4> +a* b + k utj> +iiy b ), (2.61) 

= L ^ X 4>t^K -~y 2 T v d k d4,+ )- R T k {k dl p+a*+k u<j> +iiy T ). (2.62) 

for the third family. The charged Higgs-neutral Higgs-sfermion-antisfermion couplings are 

X ^d L ui = X ^p L o* = 9 2 {K^k u ^ + fc d0 ofc d0 +)/2V2 (2.63) 
for the first two families, and 

x mht* m = L b k L U x m^i ~ ^ k *^ k u4>t +y 2 b k d^k d4) +}/V2) 

- R l k R i m ytyb{K^k d ^ + k* d4> ok u<p +)/V2, (2.64) 

X 4>1<t>p k K = L ^( A <A°0+^s* -y 2 k d 4,°k d(t> +/V2) (2.65) 

for the third family. 

The sfermion-antisfermion-sfermion-antisfermion couplings in the Lagrangian are written as 

- C = \ X hffhf^^M*rnl (2-66) 

where each combination in parentheses forms a color singlet. Then 

3 



n=l 

where the non-zero Yukawa f-term contributions are: 

f 

•v. " a 

x fir*T k f^ = VriL^R^RfkLp + RfiLf.L^R^), (2.70) 



X Wki , m = yUL u RIR ik Ll m + R~ u L\L h R\ m \ (2.68) 

m:ms. - y 2 b( L b R l R h L L + R b. L i 3 L b R V' ( 2 - 69 ) 
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X Ubfb k t* m - X b k i* m ul* - V^uH^H^+vlLhRt^H 1 *^ 
Xv T f * f fe 9; — Xf k i> T f* = y T Rf. Rf k , 

X b A ^ m = * WSt S,5j - MAL-^R^L^ + R- b LlL fk RtJ, 
X hhf k o* = X f k o*ub- = ( X b 3 t' z ^f'T = (^fjSitf )* = VbVrLt.Rl R? k - 



(2.71) 
(2.72) 
(2.73) 
(2.74) 



The electroweak J7(l)y D-term contributions to 
cq. (2.67) are: 

x' fl , = Y-j (2.75) 
for the sfermions of the first two families, and 

x' lftl = L ij LijQ-2R ij Rl/Z, (2.76) 

x' fjf . = -L fj L%j2 + R fj R* fk , (2.78) 

for the third family sfermions. The SU(2)l -D-term con- 
tributions to eq. (2.67) are 

^=^-&=e : =lA (2.79) 

*Z t = ~4U = *%l = "A- - <A (2-80) 
x (3) --s (3) -a; (3) --x (3) - 1/2 f2 81) 

for the first two family sfermions, and 
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t 3 b k 




= L h L tJ 2 ' 


(2.82) 


r (l) 
T 3 




= L%J2, 


(2.83) 


x (2) 






(2.84) 


T (2) 




= iL%j2, 


(2.85) 






= L h L V^ 


(2.86) 






= 1/2, 


(2.87) 



= -Lf 3 L%j2. (2.89) 



for the third family sfermions. The SU(3) C D-tcrm con- 
tributions to squark-antisquark-squark-antisquark cou- 
plings arc not included above, and will be kept track of 
separately in the following. 

I conclude this section with a few other important con- 
ventions to be observed throughout this paper. The sym- 
bol / refers to a generic sfermion mass eigenstate. The 
symbol q refers only to the 8 first and second family 
squarks, (u L ,d L ,u R , d R , c L , s L ,c R , s R ), which are always 
assumed to be mass eigenstates. Indices i,j are used for 
the external Higgs scalars. Indices k,m,n,p, . . . for vir- 
tual particles are always implicitly summed over all pos- 
sible values, namely 1, 2, 3, 4 for neutral Higgs scalars and 
neutralinos, or 1,2 for charged Higgs scalars, charginos, 
top squarks, bottom squarks, and tau sleptons, or over 
the 21 distinct sfermion mass eigenstates The sym- 
bol rif or refers to the number of colors, and is al- 
ways equal to 3 or 1 in the obvious way. The name 
of a particle is always used in place of its renormalized, 
tree-level squared mass when appearing as the argument 
of a loop function, so for example yL SFS p F (tk, t, t m , t, g) 
means M SFS -p F (m 2 , m 2 ,ra 2 ,m 2 ,m~). Each of the in- 
tegral functions also has an implicit dependence on s and 
Q, as in ref. [40]. All of the couplings and masses appear- 
ing below are tree-level running DR parameters in the 
MSSM with no particles decoupled. 



III. ONE-LOOP CONTRIBUTIONS TO HIGGS SCALAR BOSON SELF-ENERGIES 

In this section, I review the known results for the one-loop self-energies of the Higgs scalar bosons. The Feynman 
gauge versions of these formulas can be found in ref. [20] , but here we use the Landau gauge results in order to agree 
with the two-loop calculation of the effective potential. 

For the neutral Higgs scalar bosons 0? = (ft , H°, G°, A ), 

n i>° = \ X <t>° l <l> 3 4>l4>l k s{<Pl) + A 0»0o + - A 5 (0+) + n/^o o//.A s (/) + A o + - A o + -B 55 (</>+, </>+) 

f 

/,/' 

+2Re[Y e + d -^Yi+ 6 - (j) o]^FF(Ck,C m ) + 2Re[Y e + e -^Y e + e -^]m ek m 6m Bp T (C k ,C m ) 
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+6Re[Y m Y t ^}B FF (t, t) + 6Re[y^oF^o]m t 2 B^(i, t) + 6Re[Y b - bif)0 Yt^}B FF (b, b) 
+6Re[Y b - br Y b - b ^}m 2 b B TF (b,b) + 2Rc[r T70 or;_ o]B w (r, r) + 2Re[Y TT4> oY TT4> o]m 2 T B FF (r,r) 

+ ^9zz4> a z ^ A v( z ) + 9ww<p°4>'> A v(W) + 2Rs[g w ^+g* w ^ 4> +]B S v{4 l t ^ w ) 

+9z4,°4, ( l9z4> ( >4>'l^sv(<Pk^ z ) + ^zz^dzz^vviZ, z ) + 9ww^>l9ww<i> aB vv^-, W). (3.1) 

In general, this is a 4 x 4 matrix. It has the form of two block-diagonal 2x2 matrices in the special case of no 
CP violation. The couplings here are as defined in section II. The rcnormalizcd and finite loop-integral functions 
A s (x), A S s(x,y), B S s{x,y), B FF (x,y), Bpp(x,y), A v (x,y), B S v(x,y), and B vv (x,y) are explicitly functions of the 
tree-level squared masses of the virtual particles in the loops, and they are all also implicitly functions of s. They can 
be found in section III of ref. [40]. 

For the charged Higgs scalar bosons cj>f — (G ± ,H ± ), the result is a 2 x 2 matrix: 

f 

+ A e <t>i <t>i A e *t 07 Bss U>t > O + E n iht fi" x lt ff" Bss (/ ~' ft 

hi' 

+ ^ Y C^N m 4>+ Y 6-N m 4,+ + Y S+N m 4>r Y c+N m 4>-) B FF(C k ,N m ) 

+ ( Y C-N m <Pt Y C+N m <l>- + Y d+ N m <f>- Y 6- N m <t>^ m C k m N m B TFi C ki N m) 

+3(*atf*W+ + Y bt<t>~ Y bt<f>7 ) B FF(t, b) + 3(1^+^- + Y l^- Y tH^ mbmtE FF^ 6 ) 
+ Y r^- Y r^- B M0,r) + X -g ZZ4>u -k v {Z) + g ww<j>Uj A v (W) 
+e 2 S tJ B sv ((f>+, 0) + +B S y(0°, W) - 9 Z 4>^-9z4>+4,T B sv(<t>t ,Z) 

+9wa<pI9wa<p+ B vv(0, W) + g wz<p +g WZ(j) +B vv (Z, W). (3.2) 



IV. TWO-LOOP CONTRIBUTIONS TO NEUTRAL HIGGS SCALAR BOSON SELF-ENERGIES 

In this section, I present analytical formulas for the contributions to the two-loop self-energies of the neutral Higgs 
scalars. These are labeled in the form n^ 2 ^, where N is used to distinguish the various contributions and will be 
equal to the equation number. 



A. Strong contributions 

The contributions to the neutral Higgs scalar boson self-energy matrix involving the gluon are: 
n S*S = 4 53 2 { (2Re[F tT0 or t ^o]G FF (i, t) + 2Re[Y t - tr Y t j r ]m 2 G ¥1? (t, t) 

+X <t>°ht^<t>° j i m tl G Ss(tk,tm) + ^o^of^W sssv(tk,tk,tk,0fj +(t—>b) 

+ E X ^qq' X <k°,qq' G Ss{q-, <?) + ^ \>?tf?99* W SSSV (<?, Q, <Z, 0)}. (4.1) 
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The functions G FF (x,y), G-pp{x,y), and Gss(x,y) are defined in section V of [40]; they follow from the Feynman 
diagrams labeled M FFFFV , Vffffv (with fermion mass insertions in all possible ways) and Mssssv, Vssssv m 
figure 1 of the present paper. The contributions involving the gluino are given by: 

B % = 16 9l{[M( Y tt^L ik L* im +y t ^oi? t - fe i?| m )A^o t - mt -,]m t M SFSFF (4,t,t m ,i,3) 
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n 



n (2),4 



(2),3 _ 



-R e l Y tt<p°^ m R i k X rf>°/ k t> m ] m g M SFSFF~( t k, *> *m, *j fl) 

-R*^ji t - ra ^A o ffct % JmJmgM^^^fc, i, i ro , i, <- .7) } + (i - &), (4.2) 

oYfitjfi] [VfFFFS^i *j 9j *fc) + FF~FFs(t> t, t, g, tk)] 

+2Rc[Y m oY tJ ^]m 2 V FFFFS (t, t, t, g, t k ) - 4:Re[Y tU oY t ^o]Re[L ik R* ik ]m t mgV FFF ^ s (t, t, t, g, t k ) 
-2Rc[Y tU oY tU oL ik RlJm t mgV FFFFS (t, t, t, g, i k ) 

-2Re[F^oF t ^oi| fc i? t -Jm?m g Vyyyy s (i, t, t, g, t k )} + (t - 6), (4.3) 

8 3a{ [V^ ^^ w ssFF(t k , t k , t, 9) - 2Re[X^^o i j, m L* ik R i Jm t m s W SSFT (th,t m , t, g) 
+2Re[X^ k ^ m X^ m f,]VsssFF{tk,t m ,t m ,t, g) 

-2Re[A^o t - fc t* m x 4,°t n t* k ( L t m R i n + R t m ^J] m tm g V SSSFF (t k ,t m , t n , t, g) + (t -» 6) 

+ 51 tytfqg'WssFFiq, q, 0, ff) + 2 ^ Re[A o^, A^o^, ]V SS sff(<7, <?, g, 0, <?)}. (4.4) 



Finally, the contributions from squark-antisquark-squark-antisquark interactions proportional to g\ are: 



n 



(2),5 



4 .93( ^<P°4>° j tJ' m ( L *t k L t n - R i k R tJ( L t m L *t n ~ R t m R *t n )^SSs{tk,t m ,tn) 

+2Re[\^ jj* m \ cf) otj,(L im L^ - R im R~J(L ip L~^ - Ri p R~ n )]Yssss(tk,t m ,t n ,t p ) 
+X <t>?hi^ x <t>°t n t P ( L i m L i n ~ R t m R iJ( L t p L i k ~ R t p R i k )Zssss(tk,t m ,t n ,t p )^ + {t^b) 
+ ( x <t>°<p<>qq' x sss(q, q, q) + X^, x^, [2Y SS ss(q, q, q, q) + Zssss(q, q, q, g)]) }• 



(4.5) 



B. Yukawa and related contributions 

In this section, I present contributions to the neutral Higgs scalar boson two-loop self energy that involve Yukawa 
couplings and the corresponding soft (scalar) 3 interactions, as specified in the Introduction. 
The contributions involving charginos and neutralinos are given by: 



n 



(2),6 



= 2n t 



n (2),7 



n 



( Rc Kt0oA o t ~ mt ~, Y t *fj ki . Y iN k iJ m N k ^SFSFT ( J'rn, t, t n , t, N k 



(2),8 _ 



+Re[Y m X^ mi , n {Y*^~^ 

+Re[F^oA^o t - mt -, Y t ^ ki , m Y^ ki Jm 2 t m^M SFSFF (t m , t, t n , t, N k )) + (t <- j)] + (t - b) + (t - r) , (4.6) 
2 3^Re[y t ^oA oj m fc, y t ^ fcS . Y-* dkl Jm^M SFSFF (b m , t, b n , t, C k ) 

+ R *[ i W°^og m g. (Y^g. ^c fe ^ + Y td k i m Y tC k 0^ mtM SFSFF(bn, t, l m , t, C k ) 

+^ c \ Y tt4>« X 4,fb m b{ Y tc k l^tc + (*<-►&) 

+Rc[y rT0 o A o^,y* afc .,y_* afc .Jm 6fc M iSFSF ^(i> T , t, z> r , r, Cfe) 
+Re[y r70 oA o i>p ,](|Y T(5 , fc ^| 2 + \Y T c k0 f)m T M SFSF - F (v Tl T,v T ,T,C k ) 

+Re[Y TT4> o X^Y TdkK Y* ^]m 2 T m d M SFSFT {v Tl T, u T , r, C k )] + (i <-► j), (4.7) 

) + m 2 \ FFF : FS (t,t,t, N k ,t m )} 

+ ^4 Y tU oY tt4>'>} m t Y FFFFs(^ *' ^fe) *m) 

+2m t m^ fc Re^oF^o y t ^.t^^ fc t m ] V FFFFs(^ *> *> ^fc. f ™) 
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n 



(2); 



n 



(2), 10 



tt(2),11 



+^e[Y tJ ^oY t j^o]Re[Y t ^j^Yj^jjV FF pp S (t, t, t, N k ,t m ) 

+Rc[Y ar Y tU0 Y t ^ ki , m Y WJ jm 2 v Fms ^ t, t, N k , i m )] } + (t — 6) + (* -> r), 



(4.8) 



2n, 



(l y tc fe ^l + fed 2 )(Re[F^oy t ; o]{V FFFFS (i,i,t,C fe ,6 ro ) 



{ V FFFFS (t, t, t, Ck, b m )} 



+2Re[Y tl4> oY tJ(j> o}m t Y FFFFS {t, t, t, C k ,b m ) ) 

+2m t m 6k (M Y t^ Y tU° Y tc k h m Y tc k lJ Y FFFFs(^ *> *. &m) 
+2Re[F t ^oF t | o]Re[F i p fc ^yj ( 5 fe g m ]V FF jrp iS ,(t, i, i, C k ,b m ) 

+ Re Kt0O^ y tC fe 6^ y tC fc 6 m ]WVjWFs(*' *> *i + (* ^ & ) + (* ^ T ) . 

n t { 2Re [ A o t - m t -. A^o t - p t % (Y*^ t -. y t ^ t -, + ^ fn ^ fj> )] Vsssff (im , tn ,t P ,t,N k ) 
+ fo^*W B + Y lN k t m Y t*N k i* n )^tm^W SSFF (t m , t n , t, N k )} } + (i ^ b) + (t -+ t) 

+ \ Y vN k v* I 2 [2A o pi> , A o PtP , V SS SFf(£t, *>r, ^r, 0, Affc) + X^o P0 ,W S S F F (l>T , V T , 0, Nfe)] , 

3{2Re[A,o 5m5 , A,^ p ^ ^5.^85 + Y ic k b Y W k bJWsssF F (b m , K, b p , t, C k ) 

+ mC \-hrb m K X 4>fb p b' m ( Y tC k b n Y tC k b; + Y W k b Y hK^ mtm C* Y SSSFF(brn, K, bp, t, C k ) 
+X <Wimh [( Y tC k l. m Y ;c k ~ K + Y ic k b m Y tC k ~ bn ) W SSFF(b m , K, t, C k ) 

+ ( Y tc k b* m Y tc k b n + Y {c k b m Y tKh^ mtm c* W s^ lm > ~ bn > *' dk) ] } + {t ~ b) 
+2\ <p o j>i> ,\ (t) o j>i> ,(\Y Tdki/ ,\ 2 + \Y Tdk0T \ 2 )YsssFF{v T ,v Tl i> T ,T,C k ) 

+±\ 4> pv*^vv*^\ Y rC k O Y TC k v^ 

+X m M* [{Kd kK I' + \ Y ~C k V T I") W SS Ff(Vt, i>r, T, C k ) 

+^e[Y Tek -,Y Tdk -Jm T m e W SSFF (i> T , D t ,t, C k ) 

+2Re [ A o f mf , \$o fpf ^Y*^ ekf p' UT c kf ,]VsssFF(T m ,Tn,f p ,0,C k ) 

+ A *°^^f;^ T c t f; 5 C r c t f ,Ws5FF(r m , f n , 0, C k ). 
The contributions involving virtual Higgs scalar bosons and third-family fermions are: 



(4.9) 



(4.10) 



(4.11) 



n (2)42 



4n t 



M Y U^Vt4 ! 2 { V ^^^( £ ' *' *' *> 4>k) + m t Y FFFFS^ *> *> *> 4>l)} 



+2Re[y t ^oF^o] \Y m \ 2 m 2 V FFFFS (t, t, t, t, <jP k ) + Re[Y tir Y tI ^(Y t ^ f}m 2 Y FFFFS (t, t, t, t, eg) 
+2Re[y t ^oY t ^o]Re[(y ti0 o) 2 ]m2V FF ^ s (i, t, t, t,<f>° k ) + Re^y^ ( Y tu^ m t Y FFFFS^ *> *> *. 
+R*[***9 i t&»9 (yi ) 2 ]M FFFFS (t, t, t, t, </»°)/2 + Rc[Y tU oY*}Re[(Y t - t ,o ) 2 }m 2 M FFFFS (t, t, t, t, cjP k ) 



n ( 2 ). 13 



+Re[y^oy t ^o]|y t ^o | 2 m t 2 M FFFFS (t, t, t, t, 0g) + Re[y^oy t ^o]|y^o | 2 m 2 M FFFFS (t, t, t, t, 4P k ) 
+Rc[y t?0 oy t ^o(y t?0 oJ 2 ]m4M FFFF5 (t,i,t,i,0°)/2] + (t -> 6) + (t -» r), (4.12) 

+2Re[y ti0 oy t ^o]m 2 V FFFFS (t, t, t, 6, 0+) + (t ~ b)} 

+ 2r t60+ r 6t0- m t m 4 Re Kt0« y tt0°HV FFFFS (i, t, t, b, 4>+) + m 2 V FFFF5 (t, t, t, b, <(>+)} 
+2Re[y t ^oy^ o ]V FFFFS (i,i,t,6,0+) + (*<-> 6)} 



+21^- {Re[Y TT ^Y; T ^}[V FFFFS (T, r, r, 0, <{>+) + m 2 T V FFFFS (r, r, r, 0, <f>+)} 
+2Re[Y TT ^Y TT ^o]m 2 .V FFFFS {T, r, t, 0, </>+)} , 
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(4.13) 



n 



(2) : 14 _ 



= 6 



+ Y -bt4,- k ) m tm b {Re[Y tI ^Y* u ,}M FFFFS (b, t, b, t, <f>+) + R C [Y tI ^Y bb<p o}M FFFFS (t, b, t, b, <f>+)} 
+ Y tHl Y b^- k {MYu^Y^o] [M FFFFS (t, b, t, b, 0+) + m^mlMppppsit, b, t, b, </>+)] 

+MY m Y^][m 2 t M^ FTs (b,t,b,t,4>+) + m 2 M FFFFS (t,b,t,b,^)}}] + (i ~ j). (4.14) 
Contributions involving virtual Higgs scalar bosons and third-family sfcrmions are: 



n 



(2), 15 _ 



= n 



n 



(2), 16 _ 



^?t™t;^0°t,t* ^4> k t p t'^4,° k t n t^ M SSSSs(t m , t n , t p , t q , (f> k ) + \^o^oj mi ^\^o^oj ni , m S S Ss(4 l k> t rn,tn) 

+2Re[A^o t - mt -, X^^X^ojj, X^o^jVsssssitm, t n , t p , t q , 4>° k ) 
+2Re[(A^o t - mt -» X^o^ojj^ + A o t - mt -» A^o^ )X (j}0 jj,}XJsss S (t m , t n ,t p ,cj) k ) 

+ X 4>y° J t m t' n ^<pltJ;^4>itJt n W ssss{t m ,tn,tp, 4>k) + 2 X <py^t m t^^°4,° k tJ^sss{t m ,tn, 4>k) 
+M^j m t^4>«jA X KKU^ Y ssss(t m , t n ,t P , 4>k)] + (* -» b) + (t -» f ) + (t -» />r), (4.15) 

n t [(^?tmt; A ^6 9 S; + X <t>ft m i'\«b q b* n } X 4>tbJ* m X l+bJ; M SSSSs(t m , &n, *p, &g, 0fc ) 

+ 2Re [ A 0?t m t* X 4>°t p t* m X l+b q t* n x <t>+b q t p } Y sssss(t m , i n , t p , b g , <j>%) 
+2Re[A^og m b» X 4>fb p b* m Kt>+l p t' q X 4,+ bJ*] Y sssss(b m , b n ,b p , t q ,(j)+) 
+2Re[(X^o imi , n X 4> o (j> + ipi ^ + A^o t - mt -. A o ?i +b pt - ;k )A*+ Bpt -,]U S :sss(im, *„, &p, ^ ) 
+ 2Re [( A ^b m 6* A 0o +6 mf; + A^b m5 ,A;^+ Bmt - ; )A^^ 

+ A 0> 3 °t m t; A 0+6 p t^ A 0+fc p t* W SS55^m, ^p, 0fe ) + \,^% m ^ A^+j nf . A* + ^j, W SSSS (6 m , £>„, t p , 

+ A 0°0 3 o t m t; A 0+0-t„t^ x sss(?m, 0fe ) + 2Re[A o fmt -, A^o^^A^+^-^j^Yssss^m, tp, ^ ) 



+Aj,o j,oi_ r. A 



^b^^^-b^b 



b~ b ' x sss(L,l,^) + 2Re[A rtJ ,A,„ y A 



g.]YsSSs(bm,bn,&p,0fc) 



+2Re[A,„ 



b m t* d>°ct>ib m t* 



T,]Ssss(<f>t,b m ,t n ) + (t -> i> T ,6-> f). 



(4.16) 



Finally, the contributions involving only virtual sfcrmions are given by: 
T (2).i7 _ 



n 0o = A o o ; fc ^ [n /fc n /n A /m /; ^ j. + n~ h X~ fm - fn j, ] X sss (f k ,f m , f„ ) 

+me \- X 4>° i M* m X 4>° j U' k ( n h n h x f m kU; + n /, A / m /;/ P /;)] Yssss ^ fc ' /"» /»> A.) 
+A 0?A/^\>°/«/ p * ( n A n /» A /»/:w; + n f k X f m kf P f^ Zssss ^ k > f m > f n > 



(4.17) 



This expression includes the contributions for the sfermions of the first two families, which only have gauge interactions. 
In the numerical results of section VI, only the third- family sfermion contributions from cq. (4.17) are included. 

V. TWO-LOOP CONTRIBUTIONS TO CHARGED HIGGS SCALAR BOSON SELF-ENERGIES 

In this section, I present analytical formulas for two-loop contributions to the charged Higgs scalar boson self- 
energies, as specified in the Introduction. They are labeled in the form Tl^l'^L, where N is the equation number. 



A. Strong contributions 

The contributions to the two-loop charged Higgs scalar boson self-energy involving the gluon are: 
11$ 'l_ = 4g 2 ({Y Ib ^ + Y Ib4>+ + Y- bt<t>T Y- bt<t> -]G FF {t, b) + [Y- tb<t>+ Y- bt ^ + F^_y^ + ] mtmb G^(t, b) 
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+^4>tht* m x lfb k t* m Gss ( bk ^ m ") + X 4 > +4,jt k t' k w sssv(tkJk,tk,0) + A^+^-g fcg ,Wsssv(6 fc ,6fc,6fe,0) 

^^4>i~d L ul^'4>fd L ul [Gss(dL, Ul) + Gss(SL, Cl)] + ^-(fr+fiJqqtWsSSviQ, Q, <?) 0)) ■ (5.1) 

9 

The contributions involving the gluino are: 

= 8 53 2 A ^ fc ^{(^ -i? fmJ R? fc +Y^ t L im Ll k )mM SFgFF (b k ,b,i m ,t,g) 

+ ( Y bt4,- L t m L * bk + Y tb4>+ R i m R t k ) m b M SFSFF^rn, t, b k ,b, g) 
-( Y bt4,- L t m R i k + Y tb4,+ R i m L t k ) m 9 M SFSFFih^, t m , t, g) 

-( Y bt4,; R t m L* bk +Y^+L im RlJm b m t mgM sTs ^(b k ,b,i m ,t,g)'j + (i <-► j)*, (5.2) 
n S?J = ^{(^W+^+I^J 2 +^ r ^-I^J 2 )V™fir(t,6,6,ff,6 fc ) 
I^J 2 + *Wr I V^ffffs^ 6, 6, 5, b k ) 

+ ( Y tb4,t Y bt^- + Y bt4,- Y tb4,+ ) m tmbVFFFFs(^ b > b >9, h) 
- 2 ( Y tb4,t Y tb4,+ + Y bt4,- Y bt4,j) Rc i L h R tJ mbm 9 V FFFTs(^ 6 ' 5 > & *fc) 
- y 6t07 ^ ^ + y t fc 0+ ^b fc ^ )mf«5V 7FFfS {t,b,b,g,b k ) 

-( Y tb4>+ Y bt4>- R b k H k +Y bt<t>- Y tb^ +{t^b, 4> + <-► (5.3) 

H ?+ J = 8 ^ 2 { [V*7*V; W SSF F (t k , k, t, g) - A^-^ i? t - m + R^LtJmtmgW ssTF -{t k ,t m ,t, g) 

+ (f — > &) +^+6 fc t^^+6 fe t. [Vs5SFF(&fc,*m,im,t,3) + VssSFF^m, 6fe, 6fc, 6, ff)] 
~ X <t>t~ b ki* m K> + HK ( L im R i n + R t m L *O m t m gV S SSFT( b k,trn, t n ,t, <j) 

~\t' b ^KpJ'S L l m R bn + R l m L bJ m b m gY SSSF ^(t k ,b m ,b n ,b,g) 

+ S W SSF F {q, q, 0, 5) + x <t,+d L ui X 4>fd L ui [Ysssff(u l JlJl,0, g) 

+~VsssFF(dL, u L , ul, 0, 5) + Vsssff (cl , s L , sl , 0, <?) + Y S ssff(§l, Cl, Cl, 0, <?)]}. (5.4) 

In eq. (5.2) and in the following, the symbol (i «-> j)* means the preceding expression with i and j interchanged, and 
with complex conjugation applied to all of the couplings but not to the loop-integral functions. 
The contributions involving squark-antisquark-squark-antisquark couplings proportional to are: 



nit 5 - = ^ 



\ u - hi AL\L ln - R\R in ){L lm L\ n - R im Rl)X SSS (tk,t m ,t n ) + (t ^ b) 



+X <Ptb m ti x ly b J- k < y L l m L b p - R t m R b p )( L b n L t p - R b n R l p ) Y ssss(tk,b m ,b n ,b p ) 

+ X 4>t b kt* m X llb k hJ L t m L *t p ~ R t m R iJ( L *t n L t p ~ R i n R t p )Y S SSs(h,tm,tn,tp) 
+X 4>tbki* m Ktbnt' p ( L t k L b n ~ R t k R bJ( L i p L i m - R i p R t m ) Z SSSs(h,tm,bn,tp) 

+ X 4>+4>;qr Xsss ^ 5, q) + X 4,+d L ui x <p+d L ui [ y ssss(ul, d L ,d L , d L ) + Yssss(dL,u L ,UL,u L ) 
+Y S sss(cl,§l,sl,sl) + Y SS ss(sl,cl,cl,c l ) + Zssss(dL,UL,d L ,u L ) + Zssss(sl,cl,s l ,cl)} ■ (5.5) 

B. Yukawa and related contributions 



In this subsection, I present two-loop contributions to the charged Higgs scalar boson self-energies that involve 
Yukawa couplings and (scalar) 3 couplings, as specified in the Introduction. 



Contributions involving neutralinos and charginos are given by: 



U *t*7 = 3X ^t{^ 7 KNj t ^ 

+ ^- Y t*N n i k Y bN n b' m + Y tb^ Y tNj k Y MjJ m bM SFsTF (t k ,t, b m , b, N n ) 

+ { Y U4>- Y tNj k Y b*Nj m + Y tb^ Y t%J % Y bN n b'J m t^SFSF F (b m ,b, t k ,t, N n ) 

+ ( Y bt4,T Y bN n b- fn Y tNj k + Y tb4,+ Y t%Jl Y b\bJ m b m t m ff n M SFSFF(tk,t, l m , b, 7V„)} 

+ X 4>ff m K Y ^-KN n i>; [ Y TN n f m m N n M SFSFF^r, 0, f m , T, N n ) 

+Y T fr nf ^m T M SFSF - F (v T ,0,f m ,T, N n )] + (i <-> j)*, 
n^ 7 _ = z{{Y- tHt Y- mt \Y b ^ kl ,J + ^-^-|l5^ 5m | 2 )V FFFFS (i,&,&,^,B m ) 

I 2 + \ 2 W b v FFFFS {t, b, b, N k , hm) 

+(*W+*Wt + Y bt*7 Y tbt+WbN k l*S + \ Y w k bJ 2 )™bm t V FFFFS (t,b,b,N k ,b m ) 

+ ( Y tb^ Y b^- Y bN k b*Jm k b m + y 5^rn60+^Jv fe 6 m niv fe ^)^^iv fc V FFFFs (t, 6, 6, # fc) b m ) 

+ 2 ( Y tb4>+ Y tb4>+ + y 6t0- y 6t0-) Rc [^ fc ^ y 6JV^ m ] m t™JV fc V FFW5( < : 6 > b , N k ~b m ) 
H Y tb4>+ Y bttJ Y bN k b* m Y W k b m + Y bt4,- Y tb^+ Y lN k b m Y b\b*J m l m t m N k Y b > & > N k , 6 m )} 

+(i ~ 6, 0+ ~ 0-)* + ^ T0 r^ T 0-{l^.f m | 2 VFFFF5(O,r,r,7V fc ,f m ) 

+ l F rJV fe f^ l^r V FFFFs(°> T > T > ^> f "0 + 2Re [ y rW fc f^ y rJV fc f m ] m r™JV fc V FFFF5(°' T ' T ' ^ » f m 

+l^ T iv fc ^ | 2 Vffffs(t, 0, 0, 7V fe , £ r )}, 

n S?*7 = 3 { (y ^ F ^ +|Y ^^ 12 + F 6t C ^7l^C fet - m | 2 )VFFFFs(i,6,&,C fc ,i ro ) 

+0W+ >W+ l*fc fe * m I 2 + Y m- Y bt4>- \ Y bc k t* m \ 2 WbV FFFFS {t, b, b, c k , t m ) 

+(*W+*W 7 + *W ^^(inc.c I 2 + \ Y bC k t m \ 2 )™ b m t V FFFFS (t, b, b, C k , t m ) 
+ ( Y tb4>+ Y bt4>- Y b*c k t- m Y lc k i m + Y bt4>- Y tb4>+ Y lc k t m Y bc k h m ) m t™c k v ffffs {t,b,b,C k ,t m ) 
+ 2 ( Y tb4,+ Y tb4,+ + Y bt^- Y bt4,-)^ e \ Y bc k t^bc k tJ m b m c k ^ b > b > Ck, t m ) 

+ ( Y tb4>+ Y bt4,- Y bc k i* m Y bc k i m + Y bt<t>- Y lb<p+ Y ic k iJ Y b*c k ii) m b mtm c k ^ 

+ (t~b, Cj>+~r)* + Y T^rY T ^-{\Y Tdk ~fV F FFFs(0,T,T,C k ,i> T ) 

+ \ Y rd k o* ? m V^ ffffs^ r > r > C k ,v T ) + 2Re[F r £ fe ~, F 7(5fe - Jra r m 6 V FFFFS (0, r, r, C k , v T ) 
+ K T c k ^\ 2Y FFFFs{r,0,0,C k ,T m )Y 

+ { Y tN k hJtN k t n + nk^ y t kt*) mtm ^ W 55FF(^,i«,i,^)] + (t - 6) + (t - r) 

+ V*7 ^ N kK fWsSFF (Vt,Vt,0, N k ), 

+A tf SnC^JMS. [( y 6JV*6* + Y b*N k b n Y bN k b p ) V SSSFF(t m , b n , b p , b, N k ) 

+ ( Y bN k b* n Y bN k b p + Y lN k b Y bN k b^ m b m N k Y SSSFF&n, K, b p , b, 7V fc )]} 
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+hif nK K+f pK K^ f .^, ; + Y^ k .Y T ^ fp )VsssFF(i>r,f n , f p , r, N k ) 

+ ( Y TN k f* Y TN k f p + KN k fJrN k f ; ) m r m N k ^SSSFFiPr,f n , f p , T, N k )] 

+A^ fm -,A;+. m ., \Y VT x h -. \ 2 V S SSFF(fm,i>r, u T ,0, N k ), (5.10) 

H Y tC k b* m Y tC k b n + Y {c k b m Y ^ (* ~ 5 ) 

+\Uj*T*i [^C^A 2 + \ Y TC k ,f)W S SFF(Vr,i<r,T,C k ) 

+2Re[Y TdkP ,Y TdkP Jm T m d W ss - FW (D T ,D T ,T, C k ) 

+A^^- fmf ,y yT c fc ^ 5 C T c fc f. W ^FF(f ro ,f„,0,C fe ), (5.11) 

+ ( Y bC k i n Y bC k t p + Y bc k h Y ic k O mbm c k ^ SSSTF~( h ™, t n , tp, b, Ck)] 
+X <PtbJL X l+b P i* m tt Y tc k K Y tc k h; + Y -tc k b Y tc k b p )^sssFF{t m , k, bp, t, C k ) 
+ ( Y tc k b{ Y tc k b p + Y ic k b n Y tc k b p } mtm c k ^ sssF^fi™, k, bp, t, Ck)]} 
+ Vv;V» 5 : l(l ^ 5 f |J + l^c fcI >J 2 )VsssFF(f m ,^,^,r,C' fc ) 
+2Re[Y Tek -Y Tdk -^m T m 6 V SSS pp(f m , v T , v T , r, C fe )] 

+ht^* T X lp P ^c k f *Kc k f; VsssF F (i>r, f n , fp, 0, C k ). (5.12) 
Contributions involving virtual Higgs scalar bosons and third-family fcrmions are: 

ng£ = 3{(r Th0 +y T< 

+2(^+^7 + ^ r ^+)l^o| 2 m 6 m t V 7ji . 7ji . s (6,t,t,t ) 02) 
+2PW+ i W+ + ^t0-^t0 7 ]Re[(y t ^o) 2 ]m2V FF ^ s (6, t, t, t, 0g) 

+[*W+ 5 W 7 ( r ^) 2 + ^ r *W+ (^o) 2 ]"»6"H 3 Vtfffs(6. *- *> *> $)} + (* - fe, + - 0-)* 

|F T - o | 2 [V FFFFS (0, r, r, r, 4) + m 2 T V FFFFS (0, r, r, r, 4)] 
+2Y TVt4> - Y T ^-Rel(Y T -^f]mlV FFWF - s (0, r, r, r, 0°), (5.13) 

U *t£ = ^ Y tb,t Y tb^ t + Y M,-^ 

+t I W+ 3 W+^- + ^ t0 -^0-n 2 0+H 2 V FFFF5 (i, 6, 6, t, 0+) 

+PW+^ 7 + r 6 t 0- y ^;Hn 2 + + Y it^ m b m tV FFFFS (t, b, b, t, <f>+) 

+ + ^^- y t60+]^+ y 6t0-™ 2 [ V FFFF s( < : b > b > *> 0fc ) + r^FFFFS^ b > b > *> 0fc )] 

+2[*W+ 5 W+ +^-^7]%^+^- m 6 m t V FFFFS (i,6,6,t,<A+)} + (t <- 6, 0+ <- 0") 

F^-Fi^- [V FFFF5 (0, r, r, 0, 0+) + V FFFFS (r, 0, 0, r, <f>+)], (5.14) 

n Sj = 3 { (F ^ y ^7^oi;^o +n t0r ^;^oy^o)M™ s (t,t,6,6,^) 

+ 2 ( i a*+ i W+ + y 6t0" y 6^7 )ReK ¥ /^»] m ' m ' M ™(^ *. 6. 02) 
+2(^+%4 + Y bt*- Y bt4>- )M Y tt^ Y ^o]m b m t M FFFFS (t, t, b, b, <f>° k ) 
+ ( y th0+ *W 7 Y ii^ k Y bw° k + *W ^4 YtJ K Y bH° k ) m l M FFTFs (t,t,b,b,(f> k ) 
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HY- tHt Y- bt(j> -Y m Y; i(j>l + Y lt ^-Y- tHt YyY bl ^mlM FFT¥s {b, 6, t, t, 0°) 

+ (^0+^;^^^ (5.15) 
Contributions involving virtual Higgs scalars and third-family sfermions are: 

n i+07 = n A (V^f*. ^+07*™*; x l+bj ; v ssss(t m , *p> ^ ) + A *j0+6 p t S , X 4>lb P K V ssss(t m , K, b p , 4>l)] 

+ A 0+h m t* [ A *jf*7MS. A ^bp** Uss5s(&m, *n, &p, <Afc) + A 0» fe 0+h m t; A *°t„t; U SSSs(&m, *n, *p, <A°)]) 

+ (i <-> j)* + <V+ - t - mt % \+ <t) jt n h m ^SSs{(l>k,irn,in) + ^+0"fc m fc* A 0+07 6„6^ ^SSs{4>k , &m, M 

+ A 0^+fc m t* A 0^0+fc m t* S SSS(^°, &m, in) + A 0+fc m t; A 0+6„t- A 0°6„fc^ A 02*p*5 M SS\S\S\s(&m, & n, tp, tj, 0°) 

+X 4>t4>~t m ti [ X 4>lU; X 4> a k t p t' m w ssss(t m , t n , t p , 4>k) + x 4> +b p t* m x l+b p t* ri w ssss{t m ,t n , b p , 4%) 

+ X <t> a k <i>l in i* m X SSS (*m , *n , </>fc ) /2 + - ^ ^ X S SS (*m , *n , 0fc )] 

+ A 0+07& m fc* [ A ^^b; A 0^P^ Ws ' s ' s ' s '^ m ' 6 "' 6 P' ^) + A ^b"*; A 0^ m t; Wssss ^ m ' 6 ™' *f ^ 

+X <t>t~bni* m Ki,+b p t' m [ X 4>lb P bi X 4>ll q l* n w sssss{t m ,b n , b p , b q , <j>° k ) + A + 5pt -»A*+^ f ,V S ssss(im, &„, b p , t q , <t>t) 
+ X 4 > l4 > « k b p b fl Y ssss{t m ,b n ,b p , 4>° k )/2 + X^+^-i^Yssssitm, b n , b p , 

+A 0+fc m t* X l+h m t; [ X 4>lUl X 4>li q i; v sssss(b m ,t n ,tp, t q , <j>° k ) + \+ bqi ,\l+ bqi ,\ sssss(b m ,t n ,tp, b g , <j>l) 

+ X <P k 4>° k tJ p y 'sSSs(bm,t n ,tp,4>k) / 2 + A 0+0- 1„ t; Y SSSS (^m , ^ , , 0fe )] } + {t ^ U T ,b ^ f). (5.16) 

Finally, contributions involving only virtual sfermions are given by: 

(2) 17 ~ ~ ~ 

n +^7 = V*7/*/^^ n /* n /» A /m/:/»/; +n h X f m hfnf'J x sss(fk,f m ,fn) 

+{ - X 4>t U* m x l+ hk + x 4>t U: x l+ f m f^ n h n h x U*Jvf; + n f\ x uf;fJ^ Yssss{ J k ' f™ & 
+x <t>thf^ x l+ U; ( n h n f» x f m f:Mi + n h X f m f;krJ 7 ' ssss ^ ^ ^ (5 ' 17) 

This expression includes the contributions for the sfermions of the first two families, which only have gauge interactions. 
In the numerical results of section VI, only the third-family sfermion contributions from eq. (5.17) are included. 

VI. DISCUSSION AND NUMERICAL EXAMPLES 

I have carried out several checks on the above expressions. First, the quark/gluon two- loop contributions to the 
Higgs scalar boson self-energies had already been computed in refs. [52, 53]. I have checked that my corresponding 
results, namely the Gff and Gp-p terms in eqs. (4.1) and (5.1) of the present paper, agree with these, by converting 

from the on-shell scheme used there into the DR scheme used here. 

Second, I have verified that the self-energy contributions listed above for h°, H°, when evaluated in the limit s — > 0, 
do correspond precisely to the second derivatives with respect to v Ul Vd of the appropriate terms in the two- loop 
effective potential [31], according to: 

n( 2 Vm = I ( Ca - S A ( d2v(2) / dv l d 2 V^/dv u dv d \ ( c a s a \ 

[ ' 2 [ Sa Ca J \d*VW/dv u dv d d^/dvj ) \-s a c a ) ' ^ > 



where c a — cos a, s a — sin a, and the two-loop effective potential is 



v ^ = v m + _l_ vm + _±_ v m + ..,, (6 . 2) 



17 



Third, I have checked that the renormalization group 
scale invariance of the pole masses is consistent with the 
known two-loop renormalization group equations for the 
Lagrangian parameters [38, 54-56] and VEVs [31]. These 
checks are quite involved, but follow the same pattern as 
given explicitly in the toy model of section VI of [40] . 

Finally, there are non-realistic limits of the MSSM 
in which a global SU(2) symmetry implies the equality 
of masses and self-energies of the charged Higgs scalar 
bosons G ± ,H ± with two of the neutral scalars. This oc- 



curs for y t 
cither: 



y b , a t = a b , 



m d ., and 



case 1: 
case 2: 



v d = 0, 



9 = 9' = 0; 
9,9' 7^0, 



and neglecting all slepton contributions. I have checked 
that in each case, the required equality between neutral 
and charged Higgs scalars for the self-energy contribu- 
tions of sections IV and V indeed occurs. 

The most important application of the results above 
is probably to the calculation of the "momentum- 
dependent" contributions to the pole mass of the lightest 
scalar Higgs boson, h°. Before reporting some numerical 
examples, it seems worthwhile to illustrate the role and 
rough size of the effects with a simple limiting case that 
can be treated analytically. Consider the degenerate de- 
coupling limit in which the top squarks and the gluino 
have the same mass Af, with s <C m 2 <C Af 2 , and with 
all bottom, tau, and electroweak effects neglected. Then, 
at one loop order: 



where 



£ [Pi + *Pl 



+ . . 



Pi 
Pi 



6Af 2 [lnM 2 - 1] - 6m 2 [Era 
+121n(M 2 /m 2 ) 



1] 



31nm 2 + 2, 



and 



InX = \n(X/Q 2 ), 



(6.3) 



(6.4) 
(6.5) 



(6.6) 



and we consistently neglect terms of order m 2 /M 2 . Sim- 
ilarly, at two-loop order, we obtain from the results of 
section IV A, and the analytical expressions of section VI 
of ref. [41], and for simplicity keeping only terms of order 

9h 



where 



ni 2 oU*) = 33 2 y t 2 ^ [P2 + SPZ + ...] 

P 2 = 32Af 2 [-(hiM 2 ) 2 + 3foM 2 - 3] 
+16m 2 [9(hrm 2 ) 2 - 9Em t 2 + 5 



(6.7) 



-2(lnAf 2 ) 2 - 61nAf 2 lnm 2 + 51nM 2 ] , (6.8) 



P' 



-12(lnm 2 ) 2 



— , 44 
121nm t 2 + — 



-4(lnAf 2 ) 2 + -InAf 2 + 81nAf 2 lnm 2 . (6.9) 
3 

Now, including the tree-level contribution to the squared 
mass, one can use the condition 



dV eS /dv u = 



(6.10) 



to eliminate the terms proportional to M 2 in the expres- 
sion for the pole squared-mass. One then finds: 



"»fco, P oie = «4cos 2 (2/3) + ^c 2 [m 2 A 1 + m^A'J 

+ ^KA 2 + ^^], (6.11) 

neglecting terms of order yf and m^ /m 2 , with 

Ai = 121n(Af 2 /m 2 ), (6.12) 

A; = P{, (6.13) 
A 2 = 32[3(lnm 2 ) 2 - lnm 2 - 1 - (InAf 2 ) 2 

-2hiM 2 W 2 + hiM 2 ], (6.14) 

A' 2 = P> 2 . (6.15) 

Choosing the renormalization scale Q — M, 



Ar = 

a; = 

A 2 = 

a; = 



12L, (6.16) 

2-3L, (6.17) 

96L 2 + 32L-32, (6.18) 

-12 J L 2 + 12i + 44/3, (6.19) 



where L = ln(M 2 /m 2 ), and as usual the masses and y t 
and c/3 are DR couplings in the MSSM (with the top 
quark and the superpartners not decoupled). The terms 
Ai and A 2 agree with the results obtained in eq. (21) of 
ref. [22]. The last term, A 2 , is a consequence of the new 
result obtained under much more general circumstances 
in this paper. However, even in this crude limit (which 
neglects the important ingredients of top squark mixing 
and mass hierarchy), we can see that it is smaller than 
one might perhaps have expected. This is both because 
the dimensionless number coefficients in the A 2 term are 
smaller than those in the A 2 term, and because there is 
a significant cancellation between the leading logarithm 
squared term and the sub-leading logarithm and constant 
term in A' 2 . Indeed, the leading-logarithm approximation 
to A' 2 is clearly quite poor unless M is over 1 TeV. 

For more precise results in realistic models, it is neces- 
sary to keep all of the terms in the two-loop self-energy, 
and evaluate the integrals numerically. When computing 
the pole mass of h°, it is best to use the following trick for 
approximating the full two-loop self-energy. Denote by 
IIp ar (s) the sum of the 2x2 matrix self-energy contribu- 
tions for the neutral Higgs scalars h°, H found in section 
IV. (From here on I only apply the general results above 
to specific examples without CP violation.) Then we use 
the following expression for the two-loop self-energy: 

H( 2 )(,s) « H( 2 a ) (s) H( 2 a )(0) + H( 2 )(0), (6.20) 
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where the last term is given exactly by eq. (6.1). In 
this way, we include all other two-loop self-energy effects 
within the effective potential approximation, while avoid- 
ing any possibility of double-counting. Eventually, when 
all of the remaining diagrams are calculated, this proce- 
dure will not be necessary, of course. 

For a specific quasi-realistic numerical example, con- 
sider the model defined by the following DR parameters 
at a renormalization group scale Qq — 640 GeV: 

g' = 0.36, 3 = 0.65, 53 = 1.06, 

lft = 0.9Q, 2/6 = 0.13, y T = 0.10, (6.21) 

and, in GeV, 

Mi = 150, M 2 = 280, M 3 = 800, 
Of = -600, a b = -150, a T = -40 

and, in GeV 2 , 

TO Qi, 2 = ( 78 °) 2 > ™L = (740) 2 , m% 2 = (735) 2 , 
m 2 Li 2 = (280) 2 , m 2 ei 2 = (200) 2 , 
m% 3 = (700) 2 , m 2 U3 = (580) 2 , m 2 3 = (725) 2 , 
m| 3 = (270) 2 , m 2 3 = (195) 2 , 

m 2 ^ = -(500) 2 , m 2 Hd = (270) 2 . (6.22) 
The two-loop effective potential is then minimized by: 

v u (Qo) = 172 GeV; v d (Q ) = 17.2 GeV, (6.23) 
provided the remaining parameters are: 
H = 504.18112 GeV, b = (184.22026 GeV) 2 . (6.24) 

Figure 2 shows the two-loop contribution to the quan- 
tity RepI/jO/jO (s)] — II^o^o(O) in this model, as a function 
of s. The solid line is the total calculated in section IV of 
this paper. Various contributions to this are also shown 
separately: the part coming from diagrams involving a 
top quark loop and a gluon [the Gff and Gp-p terms in 
cq. (4.1)] are shown as the long-dashed line, the part from 
other diagrams involving top (s)quarks and gluinos are 
shown as the dot-dashed line, and all of the remaining 
contributions are lumped together as the short-dashed 
line. This shows that, at least for the subset of contribu- 
tions found in this paper, the deviation from the effective 
potential approximation comes mostly from top quark 
loops involving the strong interactions, as one might ex- 
pect. The relative proportions from different diagrams 
varies rather strongly with the choice of renormalization 
scale, but the total has only a small Q-dependence. Dia- 
grams involving only squarks contribute less to the quan- 
tity Re[n ft o/jo(s)] — H h o h o(0), because s <C m|. 

The resulting pole mass of h is shown in figure 3, as 
a function of the choice 2 of renormalization scale Q. To 



2 To avoid instabilities in the effective potential approximation to 



10 C — ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 r 



> 




FIG. 2: The two-loop contributions to Rep^o^o (s)] — 
n h o/,o(0) found in section IV, for the model described in the 
text, as a function of the momentum invariant s. 



make this graph, all of the model parameters including 
the VEVs are evolved using the two-loop renormalization 
group equations [54] from the defining scale Qq = 640 
GeV to the scale Q. The two- loop effective potential 
is then required to be minimized, determining the val- 
ues of \x and b at that scale. Using these parameters as 
inputs, the dot-dashed line shows the pole mass as calcu- 
lated in the full effective potential approximation, as in 
ref. [32]. The solid line shows the improved calculation of 
this paper, using eq. (6.20) for the momentum-dependent 
self-energy. (For comparison, the dashed line shows the 
result within a partial two-loop effective potential ap- 
proximation [21-23, 28], in which all electroweak effects 
involving g,g' are neglected in the two- loop effective po- 
tential.) We see that including the s-dependence in the 
self-energy lowers the prediction for the pole mass, by 
only about 160 MeV in this model, and nearly indepen- 
dently of the choice of renormalization scale. 

The imaginary part of the pole mass can in principle 
be used to obtain the physical decay width of h°. The 
contribution from various decay channels can be identi- 
fied by isolating the imaginary parts due to each one-loop 
and two-loop contribution to the self-energy. In fig. 4, I 
show the width corresponding to the decays h° — * bb 
and h — > bbg. (Not included are spurious imaginary 
contributions of the self-energy coming from diagrams 
with Goldstone bosons, which arise because we have not 
included all of the two-loop self-energy diagrams with 



the self-energy [31], only choices of Q leading to positive Gold- 
stone boson tree-level squared masses are shown; in this model, 
that requirement limits us to Q > 568 GeV. This includes the 
geometric mean of the top squark masses, and also the scale 
where the sum of the one-loop and two-loop corrections to m h o 
vanishes. 
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FIG. 3: The pole mass of h°, computed in various approx- 
imations, for the model described in the test, as a function 
of the renormalization scale Q. In each case, the full one- 
loop self-energy is used in the computation. The dashed line 
also includes the contributions of the two-loop self-energy in 
the effective potential approximation, neglecting electroweak 
couplings. The dot-dashed line includes the contributions of 
the full two-loop self-energy in the effective potential approx- 
imation. The solid line also includes momentum-dependent 
contributions to the self-energy, as found in section IV. 



> 
CD 



CD 

I 1 

T 



■ l-loop self-energy + 2-loop effective potential 
2-loop self-energy, expansion around tree 

■ 2-loop self-energy, iterative 



600 650 700 

Renormalization scale Q [GeV] 



750 



FIG. 4: The dependence of the h° — > bb(g) width, obtained 
from the corresponding contributions to the imaginary part 
of the pole mass, as a function of the renormalization scale 
Q, in various approximations. 



which therefore includes the effects of gluon emission and 
one-loop corrections to the h°bb vertex and the 6-quark 
propagator. The complex pole mass is obtained by iter- 
ation of eq. (1.4). In contrast, the dot-dashed line shows 
the same result, but using the method of expanding the 
self-energies about the tree-level mass, as in eqs. (1.5)- 
(1.7). The latter method has a strong Q-dependence for 
the width (although it only makes a difference of at most 
a few tens of MeV in the real part of the pole mass). 
This is because the tree-level h a mass is only close to the 
two-loop mass for renormalization scales near Q = 675 
GeV. Of course, the Higgs decay width is more accu- 
rately calculated using other methods (see e.g. [57, 58] 
and references therein). 

I have checked that comparable results obtain for a va- 
riety of other MSSM model parameters, including some 
with large tan/3. As one illustration, consider the effect 
of the top squark mixing, which is well-known to have a 
significant effect on the h° mass. Figure 5 shows the de- 
pendence of the computed pole mass on the Lagrangian 
Higgs- ik-ifj coupling parameter a t , keeping all other pa- 
rameters (except [i and b) fixed to the values given above. 
Recall from the definition of ref. [3] or [31] that the off- 
diagonal entries in the tree-level top-squark squared- mass 
matrix are v u at—[iytVd- Therefore, the top squark mixing 
angle vanishes for at = /Ltj/t/tan/3 (in this model, about 
45 GeV) . Figure 5 illustrates that the part of the h° pole 
mass coming from momentum-dependent effects in the 
two-loop self-energy is at most a few hundred MeV, and 
often much less. 

In fig. 5, the maximum h° pole mass is obtained for 
negative at, which at first sight might appear to differ 
from the results obtained in refs. [22-24, 28]. The reason 
is that different quantities are being held constant while 
varying a t . In those papers, the on-shell masses are cho- 
sen to be held constant, while in this paper the running 
parameters at the input renormalization scale are held 
constant instead. The fact that these two slices through 
parameter space give opposite results for the condition 
that maximizes the h pole mass can be immediately seen 
by comparing eq. (21) with Q 2 — m~ and eq. (27), both 
in ref. [22]. 

As a numerical study of the effectiveness of the partial 
two-loop self-energy corrections obtained in this paper, 
consider the masses of the Goldstone bosons. Because the 
self-energies are obtained by expanding the Higgs fields 
around VEVs that minimize the Landau gauge two-loop 
effective potential, the Goldstone scalars G° and G ± are 
exactly massless at two loop order. This means that the 
matrices 



non-zero s.) The dashed line shows the result coming 
entirely from the imaginary parts of one-loop bottom- 
quark diagrams, but using the (real) two-loop effective 
potential approximation in order to get the kinematics 
correct by making a reasonable approximation for the 
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each have one eigenvalue. In figure 6, I show the 
tree-level, one-loop and partial two-loop approximations 
to the Goldstone boson mass quantity rv? G / \f\rnQ\ as a 
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FIG. 5: The dependence of the computed h° pole mass on the parameter at, for the model described in the text. The left 
panel is the same approximation as the solid line of fig. 3. The right panel shows the change in the h° pole mass induced by 
including the momentum-dependent self-energy, compared to the full two-loop effective potential approximation. 
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FIG. 6: The Goldstone boson mass quantity maj \J\m 1 G \ in 
GeV, in the tree-level, one-loop, and partial two-loop approx- 
imations, for the model described in the text, as a function of 
the choice of renormalization scale Q. The G° and G lines 
are not visually distinguishable at tree-level (dashed) and one- 
loop (dot-dashed) order. The partial two-loop result for G° 
is the upper solid line and for G is the lower solid line. The 
full two-loop result for both G° and G^ 1 should be exactly 
by construction, since the fields are expanded around the 
minimum of the Landau gauge two-loop effective potential. 



function of the choice of renormalization scale Q. Here 
ttlq is defined to be the lowest eigenvalue of respectively 
the first term, the first two terms, and all three terms 
with n^ 2 ) replaced by Ilpi!-, in eqs. (6.25) and (6.26). 
Here n par is the partial two-loop approximation from 
sections IV and V. The effect of the approximation we 
have made for the two-loop self-energy is seen to be of 



FIG. 7: The two-loop contributions to the real parts of the 
self-energy functions n^o^o^) and n A o A o(s) (found in sec- 
tion IV) and T1 H + H - (s) (found in section V), for the model 
described in the text, as a function of the momentum invari- 
ant s. 



order only tens of GeV 2 for the Goldstone boson squared 
masses at s = 0, and much smaller than for the one-loop 
and tree-level approximations. 

Let us now turn to the effects of the partial two-loop 
self energy corrections found in this paper on the heavier 
Higgs scalar bosons 7J ± , H°, and A . These corrections 
are typically even smaller than for h°, both in relative and 
absolute terms, in part because they have a weaker cou- 
pling to virtual top (s)quarks, but also because there are 
non-trivial cancellations. Figure 7 shows the dependence 
of the real parts of the diagonal two-loop self-energies for 
H ± , H°, and ^4°, as a function of s. Since this model 
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FIG. 8: The computed i? ± pole mass for the model described 
in the text, in various approximations, as a function of the 
renormalization scale Q. The dashed line uses the one- loop 
effective potential minimization conditions to determine pa- 
rameters used in the one-loop self-energy. The dot-dashed 
line uses the two-loop effective potential minimizations con- 
dition, and the one-loop self-energy. The solid line uses the 
two-loop effective potential minimization conditions, and the 
partial two-loop self-energy as found in section V. 



is not far from the decoupling limit, these nearly form 
an isospin doublet, so the self-energy functions have a 
similar behavior, especially at larger s. Note that the 
A self-energy has a singular threshold at = 2m t , 
due to the effects of massless gluon exchange. The di- 
agrams of the type Vffffv and Mffffv in Figure 1 
cause threshold behavior proportional to (1 — s/imf^ 1 ^ 2 
and ln(l — s/Amf), respectively. If the pole mass were 
in the vicinity of this threshold, these singularities would 
have to be eliminated by re-summation, a topic beyond 
the scope of the present paper. In contrast, the thresh- 
old behaviors of the H° self-energy at *Js = 2m t and of 
the 7? ± self-energy at ^/s = m t + mj are continuous (but 
not differentiable). In all three cases, I have checked that 
there is a significant cancellation between the contribu- 
tions of order g\yl and those of order yf. The extent of 
this cancellation depends on the choice of renormaliza- 
tion scale. 

The resulting effect of the partial two-loop self-energies 
on the H , H°, and ^4° pole masses is rather small. Fig- 
ure 8 shows the renormalization scale dependence of the 
calculated pole mass for the charged Higgs scalars. Here, 
I do not use the trick of incorporating the effective po- 
tential results as was done for h° in eq. (6.20), since the 
effective potential approximation to the self-energy is not 
close to valid for the heavier Higgs scalar bosons. Here 
the dashed line shows the result of a purely one-loop cal- 
culation, meaning that the parameters [i, b are fixed from 
the VEVs by using the one-loop effective potential, and 
the pole mass is computed using the one-loop self-energy. 



FIG. 9: The computed pole masses of A (lower pair of lines) 
and H° (upper pair of lines), for the model described in the 
text, as a function of the renormalization scale Q. The ap- 
proximations are as in figure 8. 

The dot-dashed line uses the two-loop effective potential 
to fix /i, b, but then uses the one-loop self-energy func- 
tion to get the pole mass. This is seen to remove much 
of the renormalization group scale dependence. Using 
the two-loop self-energy contributions as found in this 
paper changes the pole mass by only a small amount, 
(and actually makes the Q-dependence slightly worse). 
The change is much smaller than the dependence on Q. 
The remaining two-loop diagrams involving electroweak 
gauge couplings and perhaps the three-loop contributions 
to electroweak symmetry breaking are therefore more im- 
portant than the diagrams calculated here for this case, 
and in particular should remove most of the remaining 
Q dependence in the calculated pole mass. However, 
the remaining theoretical error is probably already much 
smaller than future experimental uncertainties [6, 59]. 

Very similar results follow for the ^4° and H° pole 
masses. They are shown in Figure 9. The same remarks 
apply here as for H . 



VII. OUTLOOK 

In this paper, I have presented partial results for the 
two-loop self-energy functions of the Higgs scalar bosons 
in minimal supersymmetry, in the mass-independent and 
supersymmetric DR renormalization scheme. In the 
case of the lightest Higgs scalar, h°, this allows an 
improved calculation of the gauge-invariant pole mass, 
which should correspond to the kinematic mass observed 
at colliders. The size of the corrections was found in typ- 
ical cases to be of order one to a few hundred MeV. This 
is significant compared to the eventual experimental un- 
certainty to be obtained at the LHC and especially at a 
LC. 



22 



To make further progress, it will be necessary to in- 
clude the remaining two-loop self-energy corrections in- 
volving electroweak couplings. This has already been 
done in the effective potential approximation [31, 32]. 
However, it is precisely for these contributions that the 
approximation s = is not always a very good one, par- 
ticularly for diagrams in which no momentum routing can 
avoid an electroweak gauge boson. Therefore, it will cer- 
tainly be necessary to include these contributions in order 
to reduce the theoretical uncertainties to acceptable lev- 
els. It also seems clear that the leading (e.g. ylg\, yfgi, 
and y® ) three-loop contributions to the h pole mass will 
be necessary, but can be included in the effective poten- 
tial approximation. These corrections can be estimated 
in a leading-logarithm approach using the renormaliza- 
tion group, as has recently been done in ref. [34]. How- 
ever, we have seen above that the non-logarithmic pieces 
are not always small compared to the logarithmic ones. 

The size of the two-loop effects found above on the 
heavier Higgs boson masses H ± , H°, and A do not seem 
to be significant compared to the expected experimental 
uncertainties. However, I have not conducted an exhaus- 
tive search of all of parameter space, and in any case the 



marginal cost in human effort to include all of the Higgs 
scalar self-energies at two-loop order is not great, once 
the two-loop self-energy for h° is included. 

Besides calculations in the Higgs sector, it will be nec- 
essary to calculate two-loop corrections for the other su- 
perpartner masses in order to interpret the results above 
in realistic situations. This issue is particularly acute in 
the mass-independent renormalization scheme adopted 
here, since e.g. the top-quark Yukawa coupling and the 
top-squark tree-level masses are used as inputs, rather 
than the physical top-quark and top-squark masses. In 
order to make meaningful comparisons with higher-order 
calculations for the Higgs masses done in the on-shell 
schemes and to future experimental constraints or (hope- 
fully) data, the two-loop mass corrections for the top and 
bottom quark, the squarks, and the gluino, at least, will 
be needed. Fortunately, these results are definitely not 
out of reach. 
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